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Abstract 

A previously conjectured set of exact form factors of boundary exponential op- 
erators in the sinh-Gordon model is tested against numerical results from boundary 
truncated conformal space approach in boundary sine-Gordon theory related by an- 
alytic continuation to sinh-Gordon model. We find that the numerical data strongly 
support the validity of the form factors themselves; however, we also report a dis- 
crepancy in the case of diagonal matrix elements, which remains unresolved for the 
time being. 

1 Introduction 



The investigation of integrable boundary quantum field theories started with the seminal 
work of Ghoshal and Zamolodchikov [T], who set up the boundary R- matrix bootstrap, 
which makes possible the determination of the reflection matrices and provides complete 
description of the theory on the mass shell. 

For the calculation of correlation functions, matrix elements of local operators between 
asymptotic states have to be computed. In a boundary quantum field theory there are 
two types of operators, the bulk and the boundary operators, where their names indicate 
their localization point. The boundary form factor program for calculating the matrix 
elements of local boundary operators between asymptotic states was initiated in p|. The 
validity of form factor solutions was checked in the case of the boundary scaling Lee- Yang 
model by calculating the two-point function using a spectral sum and comparing it to 
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the prediction of conformal perturbation theory. In [3] the spectrum of independent form 
factor solutions in the scaling Lee- Yang model and the sinh-Gordon model was compared 
to the boundary operator content of the ultraviolet boundary conformal field theory and 
a complete agreement was found. It is also possible to compare form factors to matrix 
elements of local operators evaluated directly from the boundary quantum field theory 
in a non-perturbative framework. For periodic boundary conditions, this was developed 
in [H [5]; the extension to boundary form factors was obtained in [6] and used to verify 
the results of the form factor bootstrap in the scaling Lee- Yang model by comparison to 
boundary truncated conformal space approach. 

Further solutions of the boundary form factor axioms were constructed and their struc- 
ture was analyzed for the sinh-Gordon theory at the self-dual point in [7], and for the 
A 2 afline Toda field theory in [8]. One of the present authors constructed a solution for 
boundary exponential operators in sinh-Gordon theory [9], and the solution was checked 
by computing the conformal dimensions and vacuum expectation values of the fields in a 
cumulant expansion ordered by powers of the coupling constant, which can be compared 
to known exact results. However, the conformal dimension essentially tests only the part 
constructed out of the bulk form factors, and the expectation value was checked only to 
the lowest nontrivial order in the coupling constant. 

Our aim in this paper is to provide a detailed non-perturbative verification of the 
solution presented in [9]: using the ideas of [6] we aim to compare the form factors to 
numerically computed finite volume matrix elements. However, this cannot be performed 
directly in the sinh-Gordon theory as we have no way to construct a truncated conformal 
space in this case. Fortunately, it is easy to argue that (at least to all order of perturbation 
theory) an analytic continuation to sine-Gordon model should work, and for this model 
working truncated conformal space program was developed in [TU| ITT] . This forms the 
basis of the present work. 

The paper is structured as follows. In Section 2 the boundary form factors of exponential 
operators in the sinh-Gordon model are recalled (in doing so some unfortunate typos in 
the paper [9] are also fixed). In Section 3 these form factors are analytically continued to 
obtain breather form factors in sine-Gordon theory, and summarize briefly the necessary 
ingredients to obtain predictions finite volume matrix elements. Section 4 contains our 
numerical analysis, and we present our conclusions in section 5. 

2 Boundary form factors in the sinh-Gordon model 
2.1 Boundary sinh-Gordon model 

The sinh-Gordon theory in the bulk is defined by the Lagrangian density 

1 vn^ 
C = -(d,^) 2 - —(coshb^ - 1) 
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It can be considered as the analytic continuation of the sine-Gordon model for imaginary 
coupling. The S-matrix of the model is 
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It satisfies f(9, B) — > 1 as — >■ oo, and approaches its asymptotic value exponentially fast. 
Sinh-Gordon theory can be restricted to the negative half-line with the following action 
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which maintains integrability [I]. The corresponding reflection factor depends on two 
continuous parameters and can be written as |T3] 



R(9) 
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It can be obtained as the analytic continuation of the first breather reflection factor in 
boundary sine-Gordon model which was calculated by Ghoshal in [13J. The relation of the 
bootstrap parameters E and F to the parameters of the Lagrangian is known both from 
a semi-classical calculation |14[ [T5] and also in an exact form in the perturbed boundary 
conformal field theory framework [16J. 



2.2 Boundary form factors in sinh- Gordon theory 

Here we recall the results of [9J, but with a few typos corrected in the expression of the 
form factor polynomials in f)2.12p and f )2.16p . For a local operator 0(t) localized at the 
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boundary (located at x — 0, and parametrized by the time coordinate t) the form factors 
are defined as 

i^„(4 4 . . . , 4; U 02, ■ ■ ■ , ^e-^Kcosh^-Ecosh*;.) 

for 0i > 02 > ■ ■ ■ > 9 n > and 9 X < 9 2 < . . . < 9 m < 0, using the asymptotic in/out 
state formalism introduced in [17] . They can be extended analytically to other values of 
rapidities. With the help of the crossing relations derived in [2] all form factors can be 
expressed in terms of the elementary form factors 

out (0 1 0(0) | #i, 2 ,..., 6 n ) in = F°(9i, 9 2 ,..., 9 n ) 

The general form factor solution can be written as j2] 

Qn{yi,V2 •••,2/n) 



F n (9i, 6> 2 , . . • , 9 nJ 



UiVi Yl(yi + yj) t = \ 

Kj 



X\r{9 i )Y[f{9 i -9 j )f(9 i + 9 j ) (2.6) 



i<j 



where 

y = 2cosh9 (2.7) 

The Q n are symmetric polynomials of its variables, and the minimal one-particle boundary 
form factor is given by 

r ,m = lsinh6 U (0,B) 

(sinh 9 — i sin 7) (sinh 9 — i sin 7') ' 2 

where 
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is a convenient normalization factor. The polynomials Q n satisfy the following recursion 
relations: 



K : 



B : 



Q 2 (-y,y) = 
Q n +2(-y,y,yi, ...,y n ) = 

(y 2 - 4 cos 2 7)(y 2 - 4 cos 2 i)P n (y\y l} ... ,y n ) Q n (y u ... ,y n ) for n > 

Qi(0) = 
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icosjcosj'B n (y 1 ,...,y n )Q n (y 1 ,...,y n ) for n > (2.11) 
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where 
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with the auxiliary variable z defined as a solution of y — z + z 1 ; i.e. from (12. 7p one has 
z = e e . 

Let us define the elementary symmetric polynomials by 

n n 



Y[(x + Xi ) = J2 Xn ~ l(T l ( X ^--i X n) 
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The upper index will be omitted in the sequel, as the number of variables will always be 
clear from the context. Let us also denote 
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Let us also introduce the polynomials Pn 



(fc). 



P} k) = [k] 

= [k}detM {n \k) 



n > 1 



,x n ) i,j = l,...,n- 1 



^(fc) = [i-j + fcK^X!,^ 



which are ingredients of the bulk form factor solution |18| . 

In terms of these definitions, the form factor solution can be written as 

>(*) 



Qi fc) = e 1 e 2J B n _ 1 g^ 1 + trjpm + a n Af 



(2-15) 
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(k) 

where A n is a linear combination of products of 07 with total degree strictly less than 
n(n — l)/2, and the first term is understood with the replacement 

,r« -» «f > 

and the notation 

ei = 2cos7 , e 2 = 2cos7' 
was introduced. The A-polynomials are given by 

A ( 2 k) = 

4 fc) = W^ + e^ + ^e^)^ 

|( fc ) _ /I r r "1 2 . ri,l2/.2 . 2 . r.i v J / . a -2' kB 



AT = 4 sin — [&] oier 3 + [A;] (q + 63 + [A;]eie 2 )ar ^o" 2 + 4sin —J (2.16) 

up to 4-particle level, and it can easily be extended to higher levels using any symbolic 
algebra software. 

These form factors correspond to the field 

1 e *|»(o,*) 



fc|$(o,t) 



e 2 



which is normalized to have unite vacuum expectation value. 



3 Breather boundary form factors in the sine-Gordon 
model 

3.1 The boundary form factors of multi-5i states 

The theory is defined by the following action 
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The bulk spectrum of this theory consists of a soliton-antisoliton doublet of mass M, with 
a number of breathers (depending on the coupling (3) and their S'-matrices are known |19] . 
For the full boundary spectrum and associated reflection factors the interested reader is 
referred to [20] (and references therein). For the purposes of the present work only the 
consideration of the first breather B\ is needed. 

In this model, the exact expectation values of boundary exponential operators are 
known [SHE]: 
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which also defines the normalization of the couplings fi and [is- The coupling fi can be 
related to the mass M of the soliton as follows 122 1: 
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The analytic continuation between the sinh-Gordon and sine-Gordon theory in the bulk is 
defined by 



/8tt 

Under this substitution the sinh-Gordon particle's 5*-matrix (12. ip becomes the S-matrix of 
the first breather in sine-Gordon theory: 
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The mass of the first breather is related to the soliton mass M as 

mi = 2Msin — 
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The identification of the two models can be completed by relating the boundary parameters. 
Introducing the Ghoshal-Zamolodchikov parameters [lj 
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the reflection factor of the first breather can be written in the form 
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(3.4) 



where a denotes the boundary condition parametrized by rj and Under this identifi- 
cation the relation (13. 3ft coincides with the relation between the GZ parameters and the 
Lagrangian boundary parameters derived by Alyosha Zamolodchikov [T6l ITT] . 
Comparing this to (I2.5P one obtains the identification 
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7T 7T 

As a result, the form factors of the boundary operator 



Jk/3<t>(0,t) 



can be written in terms of the sinh-Gordon form factor solution (presented in subsection 
12. 2p as follows 
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with 



i<j 



7T 
2 



i = i£0 



7T 



(3.5) 



and 



^^(Oje^^lO) 
is the exact vacuum expectation value (13.21) . 

3.2 Finite volume form factors 

We now briefly recall the formalism developed in [6j for the description of matrix elements of 
boundary operators in finite volume, specialized for the levels that consist of first breathers. 



3.2.1 Finite volume energy levels 

Introduce the bulk phase-shifts 



S 



(0) 



js(e) 
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and their boundary counterparts 



4? (0) = e 



(3.6) 



where a denotes the boundary condition. Putting the theory on a strip of width L with 
boundary conditions a and /3, the finite volume levels can be obtained by solving the 
Bethe-Yang equations [23J: 

Q j {9 1 ,...,9 n ) = 2nl, (3.7) 
where the phases describing the wave function monodromy are given by 

Q 3 (0i, . . . , n ) = 2m!L sinh 9j + ^ {5 (9j - 6 k ) + 5 {9 j + 9 k )} + 5 {a) (0,-) + (0,) 

Mi 

Here all rapidities 0j (and accordingly all quantum numbers Ij) are taken to be positive. 
The quantum numbers can be taken ordered as I± < . . . < I n ; they must all be different due 
to the exclusion principle (S(0) = —1). The corresponding multi-particle state is denoted 
by 

\{Iu...,I n }) L 
and its energy (relative to the ground state) is given by 

n 

E h .. Jn {L) = J2 mi C0S H + °( e " ML ) 
where < 9j > is the solution of eqns. ( 13. 7p at the given volume L. The Bethe-Yang 

I J j=l,...,n 

equations gives the energy of the multi-particle states to all order in 1/L, neglecting only fi- 
nite size effects decaying exponentially with L (where \x is some finite mass scale, dependent 
on the details of the spectrum). 

3.2.2 Non-diagonal matrix elements 

For non-diagonal matrix elements, it was shown in [6J that 

|<{^...,jaiO(0)|{/i,...,/„}>L| = 



i ■ • ■ i "n, 



p(9 1 ,...,9 n )p(9[,...,9'J 



+ 0(e-^ L ) 



(3.8) 



where 



p(0i,. . . ,9 n ) = det 



9Qk{9i, ■ ■ ■ ,9 nJ 



d9 t 



kl=l ri 



(3.9) 



is the finite volume density of states, which is the Jacobi determinant of the mapping be- 
tween the space of quantum numbers and the space of rapidities given by the Bethe-Yang 
equations ( 13. 71) . In general, the phase conventions used for the exact form factors and the 
finite volume matrix elements differ, so only the absolute values can be compared. Eval- 
uating the above expansion requires analytic continuation of the form factors to complex 
values of 6 which can be accomplished using the formulae given in Appendix [A] 



3.2.3 Diagonal matrix elements 

In the diagonal case, there are disconnected contributions which can be taken care of by 
regularizing the appropriate form factor as 

F 1 _ 1 (0 n + i7T + e n ,...,0i+i7T + eiA,.-.A) = (3-10) 

•in 

n ^ n n 

TT ~ • / ] ■■■ / ] Ai...i n (#i, • • • , 6' n )e il ej 2 ...ej n + . . . 

1=1 *1=1 ln=l 

where A\...i n is a completely symmetric tensor of rank n in the indices ,i n , and the 

ellipsis denote terms that vanish when taking ej — y simultaneously. The connected matrix 
element can be defined as the 6j independent part of eqn. (I3.10p . i.e. the part which does 
not diverge whenever any of the e, is taken to zero: 

F c (0 1; 2 , 9 n ) = n\ Ai... n (e h . . . , n ) (3.11) 

where the appearance of the factor n! is simply due to the permutations of the e«. The 
formula for the diagonal matrix element reads 



({/l.. I n }\O(0)\{I 1 ...I n }) L = (3.12) 
fS 1 a , E FC ({dkheAyp ai ...a n (0i, • • • , 9 n \A) + 0(e^ L ) 

PWl, AC {l,2,...n} 

The summation runs over all subsets A of {1, 2, ... n}. For any such subset the appropriate 
sub-determinant 

p(6 1 ,...,6 n \A) = detJ A (6 1 ,...,6 n ) 
of the n x n Bethe-Yang Jacobi matrix 

J(ei,...Ah = dQk{9l --> 6n) (3.i3) 

Ot>i 

can be obtained by deleting the rows and columns corresponding to the subset of indices 
A. The determinant of the empty sub-matrix (i.e. when A = {1, 2, . . . n}) is defined to 
equal 1 by convention. 

Note that in contrast to (13. 8p in (13 . 12[) it is not necessary to take the absolute value, 
as phase redefinitions drop out from a diagonal matrix element. 

Eqns. (I3.10p and (I3.12p are expected to give the finite volume matrix elements to all 
order in 1/L, neglecting only finite size effects decaying exponentially with L jU [5j |6]. 



10 



4 Numerical verification 



The finite volume energy levels and matrix elements can be evaluated using the boundary 
truncated conformal space approach (BTCSA) |2l] which is an extension of the method 
developed by Yurov and Zamolodchikov [25] to boundary quantum field theories. For the 
sine-Gordon model the program used here was developed in [TOlIll] . to which the interested 
reader is referred for details. 

To perform a specific check of the boundary form factor solution (13. 5p we choose the 
operator (corresponding to k = 1) 

Q _ e W(0,t) 

and evaluate its matrix elements between the BTCSA eigenvectors numerically. As in [6], 
all energies and matrix elements are measured in units of (appropriate powers of) the char- 
acteristic mass scale (here given by the soliton mass M), and we also use the dimensionless 
volume variable I = ML. States can be identified by matching them with the energy levels 
predicted by the Bethe-Yang equations, and the appropriate matrix elements can then be 
compared to the predictions of (13. 8 j) and A3. 12)) . obtained by substituting the exact form 
factor solution ( 12. 6p . The typical BTCSA cutoff value was 16, resulting in a truncated 
Hilbert space with several thousand states. Only a small, but representative sample of our 
numerical results are presented; we took care to verify our results for numerous different 
values of the model parameters £, rj, 



4.1 Ground state energy and vacuum expectation value 

Before embarking on the evaluation of the form factors, the accuracy of BTCSA data can 
be tested by extracting the bulk and boundary vacuum energy constants. For each bulk 
coupling £ we considered a number of different boundary conditions. In all cases, on one 
side of the strip (x = L) the boundary condition /3 is a pure Neumann one 

Vo = Vn = ~~2^^ ^o = 

and varied the boundary condition a on the other side x = 0; from now on we only give 
the parameters rj, i? of this other boundary condition a which is also where the boundary 
operator 0\ is located for the matrix element calculations. As discussed in [TO], the energy 
of the ground state is predicted to be 



E (L) = BM 2 L + (B b ( VN , 0) + B b ( V , #)) M + O (e^ L ) 



where 



1 < 
B = --tany 

Bb(v,$) = t ( cos (£t?) + cosh (£■&) cos ( — ] H — sin ( — ] 

W ' 7 2cosf V 2 V 2 J 2 V 2 J 2 / 
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t 


V/Vn 


■& 


B (exact) 


B (BTCSA) 


Bb (exact) 


B b (BTCSA) 


50/391 


0.5 





-0.050904 


-0.05009 


-0.33307 


-0.33238 


50/311 


0.7 





-0.064512 


-0.06350 


-0.16625 


-0.16802 


50/239 


0.9 





-0.085246 


-0.08355 


0.04476 


0.03796 


50/311 


0.5 


0.2 


-0.064512 


-0.06030 


-0.33289 


-0.33204 


50/391 


0.7 


0.5 


-0.050904 


-0.05039 


-0.17695 


-0.17767 



Table 4.1: Bulk and boundary energy constants from BTCSA compared to the exact 
predictions, with Bf, = Bb(r] N , 0) + Bb(r},"d). 

are the bulk and boundary energy constants in units of the soliton mass M. The comparison 
between the BTCSA results and the above exact predictions is illustrated in table 14.11 

It is also possible to test the exact vacuum expectation value against the BTCSA. The 
finite volume expectation value is expected to behave as 

(0|e w(0 '* ) |0) L = ^i + O(e^ L ) 

where the exact vacuum expectation value 

g 1 = (0|e w(0 '' } |0) 

can be evaluated from eqn. (13. 2\i . This is illustrated in figure 14.1} the deviations at 
large volume (/ > 5) are due to truncation effects. Note that agreement is better for 
smaller £, which agrees with the general trend observed from all the spectral data (masses, 
energy levels, bulk and boundary energy constants) that (B)TCSA converges better if the 
perturbing operator is more relevant (i.e. its conformal dimension is smaller )0. It is also 
important to notice that the derivative of the boundary energy with respect the coupling 
constant /i# is the expectation value of the boundary perturbation jTT], and so the boundary 
parameter dependence of these expectation values have already been implicitly checked by 
the data in table 14.11 

4.2 Level identification 

In order to evaluate matrix elements it is necessary to identify the finite volume levels. As 
illustrated in figure 14. 2\ this is performed by matching predictions from the Bethe-Yang 
equations (13. 7\i to the BTCSA spectrum. It is apparent that this is not an easy task, 
in fact, a magnitude harder than in the case of bulk theories. In contrast to the case of 
periodic boundary conditions [26], there is no conserved momentum, and so the Hilbert 
space cannot be split into sectors on the basis of momentum. As a result, the continuum 
(in infinite volume limit) starts at the one-particle threshold, as opposed to the bulk case. 
Since topological charge is not conserved either, another way of reducing dimensionality is 

1 It also happens to be the case that for a given value of £ convergence is better for smaller values of 
the boundary parameter r\. 
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^0.8 
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! T ! ! ! ! f 1 S 



. 50/391; 0.5; 0, BTCSA 

50/391; 0.5; 0, theory 

0.6 . 2/7; 0.7; 0, BTCSA 

— -2/7; 0.7; 0, theory 
a 50/391; 0.7; 0.5, BTCSA 
0-4 — 50/391; 0.7; 0.5, theory 



0.2 



5 10 15 20 

/ = ML 



Figure 4.1: Vacuum expectation values: a comparison between the exact prediction and 
BTCSA. The model parameters are listed as £; rj/rj^; All values are in units of the 
soliton mass M. The straight lines are the infinite volume vacuum expectation values from 
eqn. (13. 2p . the discrete dots are the BTCSA data. 
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5 10 15 20 



/ 

Figure 4.2: Level identification for £ = 50/391, rj = Q.brjN, d — 0. The dots are the BTCSA 
energy levels (with the vacuum subtracted), the continuous lines are Bethe-Yang predic- 
tions for three one-particle levels |{1}), |{2}), |{3}) and two two-particle levels |{1, 2}) and 
{1,3}}. 



also lost, resulting in an extremely dense spectrum. Additional complexity arises from the 
complicated particle spectrum and the presence of boundary bound states [271 120] (in fact, 
the lowest excited level in figure H~2l is precisely such a state and the next two are also clearly 
visible in the spectrum). As a result, in contrast to the bulk theories investigated in [HE] 
there is no way of reliably identifying states with more than two particles. It is also apparent 
from the figure that the two-particle states are already in a very dense part of the spectrum, 
and their identification is made harder by the numerous level crossings characteristic of 
integrable models. At certain values of the volume L there can be more than one BTCSA 
candidates for a given Bethe-Yang solution; identification can be completed by selecting 
the candidate on the basis of one of the form factor measurements, which still leaves other 
matrix elements involving the state as cross-checks. 

Level crossings also present a problem in numerical stability, since in their vicinity the 
state of interest is nearly degenerate to another one. Since the truncation effect can be 
considered as an additional perturbing operator, the level crossings are eventually lifted. 
However, such a near degeneracy greatly magnifies truncation effects on the eigenvectors 
and therefore the matrix elements |6J. 
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4.3 The one-particle form factor 

One of the most important ingredients of the boundary form factor bootstrap is the minimal 
boundary form factor ( 12. 8ft which can be obtained via measuring the one-particle matrix 
elements 

Fi(0) = (Ole^^lBi^)) 
= GxH ir (6) 

(note that the polynomial Q\{y) = P\{y) = [k]y — y for k — 1). Eqn. (13.81) can be used to 
express 

Fx(9)\ = ^/p~A^)\(0\O(0)\{I}} L \+O(e-^) 
where 9 solves the Bethe-Yang equation 

Qi(0) = 2miLsinh 9 + 5 {a) (6) + 5 {p) (6) = 2ttI (4.1) 

and 

Pi(0) = Q'M (4.2) 

Therefore it is possible to compare the values extracted from several different one-particle 
lines (distinguished by the value of J) on the same plot. The numerics and the theoretical 
prediction are in excellent agreement, as shown in figure 14.31 Note that the numerics 
deviate from the prediction for small 9 (large L) due to truncation errors, while for large 
9 (small L) the exponential corrections show up. The advantage of lower / states is that 
the scaling regime corresponds to the low-0 domain, while higher I states are useful in 
scanning the large-# behaviour of the form factor function. 



4.4 Two-particle form factors 

For the two-particle form factor, three independent tests can be performed: 
1. Vacuum-two-particle matrix element: 

^'l 1 V ■ - 17 2 . 



|(0|e i ^°'*)|{/ 1 ,J 2 }) J 



0(e 



-tiL\ 



where 6 I 1 ,6 I 2 are solutions of the two-particle Bethe-Yang equations, p n is the corre- 
sponding Bethe-Yang Jacobian, F n is given by the n = 2 case of ( 13. 5p . and the exact 
Q is inserted for proper normalization of the operator (figure I4~4"|) . 



2. Non-diagonal one-particle-one-particle matrix element: 

F n (i-K + 9',9) 

W(o,t)i- 



({/'}|e W(0 ^|{/}> 



Pi{9')pi{9) 



+ 0(e 



-liL\ 



15 




Figure 4.3: One-particle form factors extracted from BTCSA and compared to the boot- 
strap prediction. Continuous lines are the bootstrap predictions, while the circles, squares 
and triangles show data extracted using I = 1,2,3 one-particle lines respectively. The 
parameters (Cv/Vn,^) for the three figures are: (50/391,0.5,0), (50/311,0.5,0.2) and 
(50/391,0.7,0.5). lfi 
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Figure 4.4: Vacuum-two-particle matrix elements for two different two-particle states, at 
£ = 50/311,77 = 0.7777V, = 0. 



where 9 and 0' are solutions of the one-particle Bethe-Yang equations (14. ip with 
quantum numbers / and I' and p 1 is the one-particle Bethe-Yang Jacobian (14. 2 p 
(figure S3]). 

3. Diagonal one-particle-one-particle matrix element (/ = I' case): 

{{I}\e#«M\{I}) L = MlLtMj + & + o(e"^) 

(figure 14. 6p . 

The first two tests show excellent agreement!!, while the third one reveals a striking discrep- 
ancy. Since this disagreement is only seen in diagonal matrix elements, we think that the 
exact form factors are correct, and the issue is with the finite size corrections in diagonal 
matrix elements. A detailed discussion is given in the conclusions. 

2 Even so, some points are slightly displaced. As discussed in subsection l4.2l this is due to the occurrence 
of line crossings at the particular value of the volume, where identification of the state becomes more 
difficult and BTCSA cutoff errors are also magnified. 
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Figure 4.5: Non-diagonal one-particle-one-particle matrix elements for three different 
choices of the one-particle states, at £ = 50/391,7] = 0.7rj N ,$ = 0.5. 
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Figure 4.6: Diagonal one-particle-one-particle matrix elements for three different choices 
of the one-particle state, at £ = 50/391, 77 = 0.7rj N ,-d = 0.5. 
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Figure 4.7: One-particle-two-particle matrix elements for £ = 50/391,77 = O-^Vn^ — 0. 
4.5 Higher form factors 

The identified one- and two-particle states can be used to obtain tests of the three- and 
four-particle form factors from (|3.5j) . according to the formulae 



5 W(o,t) 



|{Jl,J2»L| 



+ 0(e-» L ) 



and 



({I[A}\eWW\{ Iljl2 }) L \ 



Pi(9>)p u (6 1 ,6 2 ) 

F ini (in + e' 2 ,tn + e' 1 ,e 1 ,e 2/ 



+ 0(e 



, Pii(0[,e' 2 ) Pll (e l ,e 2 ) 

(the latter is only valid for the non-diagonal case). These tests, illustrated by figures I4T71 and 
I4.8[ a l so show a good agreement between numerical results and theoretical expectations. 
Using f|3.12p diagonal two-particle-two-particle matrix elements can also be constructed, 
and they reveal the same disagreement as the one-particle ones in the previous subsection. 



5 Conclusions 

We have performed an extensive test of the form factors conjectured in |9j by comparing 
them to results obtained using the boundary truncated conformal space approach. For gen- 
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eral matrix elements an excellent agreement was found. However, for the case of diagonal 
matrix elements theoretical expectations and numerical results clearly disagree. 

The first point to make that the conjectured form factors are probably correct, as 
they survived all the tests with the exception of diagonal matrix elements. Our methods 
probed them up to four particles, at which level all of the ingredients of the form factor 
bootstrap (minimal form factors, recursion relations from bulk and boundary poles) are 
already heavily involved. Therefore we do not expect the problem to be related to the 
bootstrap. 

The finite size description, on the other hand, have only been obtained for diagonal 
scattering theories and the formula A3. 12[) for the diagonal matrix elements is only an 
educated guess based on the bulk case. One issue is that sine-Gordon model is a theory 
with a non-diagonal scattering; however, the breather scattering is diagonal, and therefore 
it is unlikely that this is the cause of the discrepancy. The other issue is the validity of the 
conjectured form of disconnected terms for the diagonal case; however, this was thoroughly 
tested both in the bulk and the boundary cases [5j [6] and there is no reason to expect any 
modification for sine-Gordon model. 

The most likely reason for the disagreement is the presence of so-called /i-terms which 
can be attributed to the fact that the first breather is a soliton-antisoliton bound states. 
Such effects were observed in [HE] and their detailed description was given in |28j . It is a 
very interesting fact that similar effects were found for bulk breather form factors (but not 
for solitonic ones!) in sine-Gordon theory with periodic boundary conditions [29], which 
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strengthens our suspicions about their origin. However, detailed investigation of this issue 
requires a knowledge of solitonic form factors of the exponential operators, which have not 
yet been constructed. We intend to investigate this issue further in the near future. 
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A Analytic continuation of minimal form factors 

The integral representations for the bulk (12.21) and boundary (I2.8P minimal form factors 
only converge in a suitable strip around the real axis. To evaluate them at complex values 
with imaginary parts as large as tt, as required in (13. 81) and ( 13. 12ft they must be continued 
analytically by separating pole factors that fall inside the strip of interest. This can be 
accomplished using the identity 



V(0;a) 



exp 

N 



dt 



t \ 2 sinh | 



cos (— ) sinh at 



sinh t 



n 



7r 2 (2fc - a) 2 + 9 2 
L \n 2 (2k + a) 2 + 9 2 

00 dt 



x exp 



t \ 2 sinh 



T - (N + 1 - Ne~ 2t ) e" 



-2Nt 



cos (—J sinh at 



sinh t 



where the natural number N is a regulatory parameter such that the value of the functions 
are independent of N, but the width of the convergence strip of the integral part grows 
with increasing N. 

The bulk minimal form factor can be expressed as 



f(0) = V(0 - m\ B/2 - 1)V(0 - in; -B/2)V{6 - in; 1) 
and for the boundary minimal form factor function u{9) (12. 9p one can use 

u(6) = U{6; B/4)U(0; 1/2 - B/A)U{6] 1/2) 

where 



;i.2) 



U(e-a) = V(9;a)V(9 -in; a 
r0 ° dt 



exp 



a 



2 cosh | cos ((z| — 9) |) sinh at 



t \ sinh | 



sinh t 
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